import numpy as np
from matplotlib import pyplot as plt

x = np.linspace(0,1000,1000)               #0表示积分下限，1表示积分上限
#y = 3*x**2                              #y表示概率密度函数
y=(1/1000)*np.exp(-x/1000)              #指数分布的例子，注意计算的是0-1000的
dx = x[1] - x[0]                        #每个矩形的宽度
fArea = np.sum(y*dx)                    #矩形宽*高，再求和
print("Integral area:",fArea)


plt.axis([np.min(x),np.max(x),0,np.max(y)])           #坐标范围
plt.plot(x,y)            #画曲线，带图示
plt.fill_between(x,y1=y,y2=0,where=(x>=0)&(x<=1000),   #填充积分区域，0表示积分下限，1000表示积分上限
                 facecolor='blue',alpha=0.2)
plt.legend()                                          #显示图示
plt.show()